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Abstract. Central compact objects (CCOs) are neutron stars that are found near the center 
of supernova remnants, and their association with supernova remnants indicates these neutron 
stars are young (< 10 4 yr). Here we review the observational properties of CCOs and discuss 
implications, especially their inferred magnetic fields. X-ray timing and spectral measurements 
suggest CCOs have relatively weak surface magnetic fields (~ 10 10 — 10 11 G). We argue that, 
rather than being created with intrinsically weak fields, CCOs are born with strong fields and 
we are only seeing a weak surface field that is transitory and evolving. This could imply that 
CCOs are one manifestation in a unified picture of neutron stars. 
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1. Introduction 

There are a variety of manifestations/classifications of neutron stars. Here we discuss 
the class dubbed central compact objects (CCOs). CCOs are very loosely defined but 
are generally characterized by the following observed properties: (1) CCOs are associated 
with supernova remnants (SNRs) and are therefore young (with ages < a few x 10 4 yr), 
(2) CCOs possess thermal X-ray flux that is relatively constant (with X-ray luminosity 
Lx ~ 10 33 ergs s~ x and a spectrum that can be fit by blackbodies from small hot emitting 
areas), and (3) CCOs have no optical or radio counterpart or pulsar wind nebula (see 
De Luca 2008; Gotthelf & Halpern 2008, for observational review, including other CCOs 
not discussed here; see also Halpern & Gotthelf 2010). 

Only three CCOs currently have a spin period P measured (as well as a measurement or 
upper limit on the time derivative of spin period P): PSR J0821— 4300 in SNR Puppis A 
hasP = 0.112 sandP < 3.5 xl0~ 16 s s" 1 (Gotthelf et al. 2010), IE 1207.4-5209 in SNR 
PKS 1209—51/52 (also known as G296. 5+10.0) has two comparable timing solutions with 
P = 0.424 s and P = 2.13 x 10~ 17 s s" 1 or 1.26 x 10~ 16 s s" 1 (Halpern & Gotthelf 2011), 
and PSR J1852+0040 in SNR Kesteven 79 has P = 0.105 s and P = 8.68 x 10" 18 s s" 1 
(Halpern & Gotthelf 2010). We hereafter refer to these three CCOs as Puppis A, IE 1207, 
and Kes 79, respectively, and only discuss them since we are primarily interested in their 
magnetic fields. 

The spin period derivative values for CCOs are low compared to most radio pulsars 
(see Fig. [l}. Assuming their current P is a historical maximum or constant (which may 
not necessarily be true; see, e.g., Muslimov & Page 1996; Geppert et al. 1999; Ho 2011; Ho 
& Andersson 2012; Pons et al. 2012), then (1) their current spin period is approximately 
their spin period at birth, (2) their characteristic age r c (— P/2P) ^ true age, where the 
true age of Puppis A is 4450 ± 750 yr (Becker et al. 2012), IE 1207 is 7 kyr with a factor 
of 3 uncertainty (Roger et al. 1988), and Kes 79 is 5.4-7.5 kyr (Sun et al. 2004), (3) their 
X-ray luminosity cannot be powered by rotational energy loss since Lx > E — 4n 2 IP/P 3 
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Figure 1. Pulsar spin period P versus spin period time derivative P. Small open circles are ob- 
served values taken from the ATNF Pulsar Catalogue (Manchester et al. 2005) . Diagonal lines in- 
dicate constant characteristic age (= P/2P) and inferred mag netic field [= 3.2 xlO 19 G [PP) 1/2 ]. 
Stars (or arrow for upper limit on P) denote CCOs at their observed values (IE 1207 has two 
comparable solutions for P). Large closed circles indicate where these CCOs would be at age 
= 10 4 yr if their initial, crust-confined magnetic field (given by the left number) is submerged 
by an accreted mass AM/Mq —Ax 10~ 4 ; right number is the surface magnetic field at 10 4 yr. 

(which could explain their non-detection in radio; see, e.g., Ho 2012), and (4) they possess 
weak magnetic fields B ~ 10 10 — 10 11 G (see next). 

2. Magnetic field of CCOs 

There are two primary methods for determining neutron star magnetic fields. The 
first involves timing observations, i.e., measuring neutron star spin period P and period 
derivative P. Assuming that the rotational energy of the pulsar decreases as a result of 
the emission of magnetic dipole radiation, the field at the magnetic equator B e can be 
inferred from P and P, i.e., B c = 3.2 x 10 19 G (PP) 1 / 2 ; note that numerical calculations 
of pulsar magnetospheres yields B e « 2.6 x 10 19 G [PP/ (1 + sin 2 a)] 1 / 2 , where a is the 
angle between the rotation and magnetic axes (Spitkovsky 2006). Using their measured 
values of P and P, CCOs have an inferred magnetic field B ~ 10 10 — 10 11 G (see Fig. [TJ. 

The second method involves spectral measurements, i.e., identifying features in the 
neutron star spectrum with particular magnetic processes. Puppis A has a possible emis- 
sion line at 0.7-0.8 keV (Gotthelf & Halpern 2009; De Luca et al. 2012), and IE 1207 has 
broad absorption lines at 0.7 and 1.4 keV (Sanwal et al. 2002; Bignami et al. 2003). If we 
assume that a spectral line at energy E is due to electron cyclotron resonance, then the 
magnetic field is B = 10 11 G (£7/1.16 keV)(l + %), where 1 + z s = (1 - 2GM/c 2 R)- 1 l 2 
is the gravitational redshift for a neutron star of mass M and radius R. The observed 
lines suggest that CCOs have B ~ (7 — 9) x 10 10 G, in agreement with the fields inferred 
from timing measurements (spectral fits of another CCO, the 330-yr-old neutron star in 
SNR Cassiopeia A, suggest it has B < 10 11 G; Ho & Heinke 2009). 
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The magnetic field of CCOs are in contrast to the majority of neutron stars, which 
possess B ps 10 12 — 10 13 G, as can been seen from Fig. [I] Furthermore, from population 
synthesis studies, neutron star magnetic fields at birth follow a lognormal distribution 
with an average and width a of logB = 12.95 ± 0.55 in the case of no field decay 
(Faucher-Giguere & Kaspi 2006) and log-B = 13.25 ± 0.6 when accounting for (model- 
dependent) field decay (Popov et al. 2010). The natural question is then one of creation 
versus evolution: Are CCOs born with weak fields or are CCOs born with strong fields 
but evolve in such a way that they appear to have weak fields at an age of < 10 4 yr? 

Halpern et al. (2007) argue for the former and propose that CCOs are neutron stars 
that are born spinning slowly. Because of their slow rotation, the dynamo mechanism 
for magnetic field generation is ineffective, and as a result, CCOs possess weak fields 
(< 10 11 G). However, there appears to be several problems with this creation scenario. 
First, CCOs do not spin particularly slowly, as illustrated in Fig. [1] this is supported 
by population synthesis work, which yields a normal distribution for neutron star spin 
periods at birth with an average and width of P — 0.30±0.15 s in the case of no field decay 
(Faucher-Giguere & Kaspi 2006) and P = 0.25 ±0.1 s when accounting for field decay 
(Popov et al. 2010). Second, a birth field < 10 11 G would require CCOs to be > 4cr from 
the peak of the neutron star distribution, and therefore there should be very few of them 
relative to the normal pulsar population. But this is counter to their observed numbers. 
For example, De Luca (2008) finds six CCOs, compared to fourteen radio pulsars, in 
all known SNRs within 5 kpc (see also Halpern & Gotthelf 2010), and Kaspi (2010) 
estimates a CCO birthrate of ~ 0.0004 yr -1 (since all known CCOs are < 7 kyr old) and 
> 10 6 CCOs in the Galaxy (comparable to the number of strong magnetic field neutron 
stars). We consider now an alternative to the creation scenario, namely evolution. 

3. Modeling magnetic field evolution 

In the evolution scenario, CCOs are born with strong fields (B > 10 12 G), and 
these fields either decayed rapidly to their current strengths or were buried by an early 
episode of accretion and are emerging or emerged recently. Magnetic field diffusion 
and decay conventionally occurs on the Ohmic timescale rohm = 47rer c L 2 /c 2 ~ 4 x 
10 5 yr (cr c /10 24 s _1 )(L/l km) 2 , where <r c is the electrical conductivity, L is the length- 
scale over which decay occurs, and 1 km is the approximate size of the stellar crust (see 
Goldreich & Reisenegger 1992). Thus fast decay from ~ 10 13 G to ~ 10 11 G could only 
have occurred in CCOs if the field is confined to very shallow layers in the star (Ho 2011). 

In Ho (2011), we compare the observed properties of Puppis A, IE 1207, and Kes 79 
to our calculations of the evolution of a buried magnetic field. We assume the field is 
buried deep beneath the surface (Romani 1990), perhaps by a post-supernova episode of 
hypercritical accretion (Chevalier 1989; Geppert et al. 1999; Bernal et al. 2010). These 
fields then diffuse to the surface on the timescale of 10 2 — 10 4 yr, so that only now 
do we see a surface field ~ 10 10 — 10 11 G. We solve the induction equation, dB/dt = 
— V x [(c 2 /47T(7 c )V x B] ~ B/rohm, in one spatial dimension (see also Muslimov & 
Page 1995; Geppert et al. 1999, for non-CCOs), while Vigano & Pons (2012) perform 
two-dimensional simulations (and thus are able to account for Hall drift) of burial and 
emergence of magnetic fields in CCOs. Fig. [1] shows examples of how evolution of an 
initially submerged magnetic field could change P and P for the CCOs. For a field that 
is confined to the crust, we also find a unique relationship between accreted mass AM 
and birth magnetic field, with a minimum logi? w 11.4 — 11.7. We find that measuring 
dB/dt or the pulsar braking index would allow a determination of AM, B, and the field 
configuration, e.g., the field is purely in the crust if dB/dt < 0, while AM is large and the 
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Figure 2. Fully ionized hydrogen atmosphere model spectra for effective temperatures 
logTeff = 5.5 - 6.7 and magnetic fields B = 10 10 ,4 x 10 10 ,7 x 10 10 ,10 n G, where the field 
is oriented parallel to the surface normal. 

field is buried deep if dB/dt > and large. We note that the (candidate) emission line 
seen in Puppis A has decreased in energy by 10% in 8.5 yr (De Luca et al. 2012); if this is 
associated with a decaying magnetic field, then a purely crustal field is implied, although 
the decay may be too rapid. We also note that optical/IR observations of IE 1207 place 
a limit of AM < 1O _6 M0 on the initial mass of a debris disk (De Luca et al. 2011). 

4. Modeling the magnetized atmosphere spectrum of CCOs 

In addition to advancements in understanding CCO timing properties, progress has 
been made in modeling their spectra. The observed thermal radiation originates in a 
thin atmospheric layer (with scale height ~ 1 cm) that covers the stellar surface. The 
properties of the atmosphere, such as the magnetic field, chemical composition, and 
radiative opacities, directly determine the characteristics of the observed spectrum (see, 
e.g., Zavlin 2009, for review). Very importantly, magnetic fields B > e 3 mlc/h 3 — 2.35 x 
10 9 G significantly increase the binding energy of atoms, molecules, and other bound 
states, and their abundances can be appreciable in the atmospheres of neutron stars 
(see Lai 2001, for review). Furthermore, when B ~ 10 11 G (T/a few x 10 6 K), models 
of atmosphere spectra must properly account for quantum and thermal effects in the 
Gaunt factor or Coulomb logarithm, which give rise to strong cyclotron harmonics in 
the opacity (Pavlov & Panov 1976; Pavlov et al. 1980; Potekhin 2010; Sulcimanov et al. 
2010, 2012). These effects are needed in order to interpret the strong absorption lines 
seen in IE 1207 as the result of electron cyclotron resonance in an atmosphere with 
B ~ 7 x 10 10 G, where the observed 0.7 and 1.4 keV lines are the fundamental and first 
harmonic, respectively (see, e.g., Ho & Mori 2008, for alternative interpretations). 

We construct fully ionized hydrogen atmosphere models using the method described 
in Ho & Lai (2001) and using Potekhin & Chabrier (2003) to calculate Gaunt factors 
and Suleimanov et al. (2012) to account for thermal effects. Examples of the resulting 
spectra are shown in Fig. [5] We note that these weak field (B = 10 10 — 10 11 G) neutron 
star atmosphere spectra will be implemented in XSPEC under NSMAX (Ho et al. 2008), 
while partially ionized hydrogen models will be the subject of future work. The spectra 
shown in Fig. [2] only describe emission from either a local patch of the stellar surface 
with a particular effective temperature and magnetic field or a star with a uniform 
temperature and radial magnetic field of uniform strength. By taking into account surface 
magnetic field and temperature distributions, we can construct more physical models of 
emission from neutron stars (see Ho 2007, for details). As an illustration, Fig. [3] shows 
the phase-resolved model spectra, pulse profile, and pulse fraction. We assume here that 
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Figure 3. Top: Atmosphere model spectra (raw 
and convolved with XMM-Newton EPIC-pn en- 
ergy resolution) at different rotation phases. In- 
set: Energy-integrated (0.01-10 keV) light curve as 
a function of rotation phase. Bottom: Pulse frac- 
tion as a function of energy, where pulse fraction 
= (-Fmax - i ? min)/(i ? max + F m i a ) . The dashed horizon- 
tal lines indicate the observed pulse fractions over the 
given energy range for IE 1207 (see De Luca et al. 
2004). 

1 + z g = 1.235 and angles between rotation and magnetic axes and between rotation 
axis and observer are 5° and 25°. The hot spot covers magnetic colatitudes — 30° and 
has effective temperature T c ft = 2 x 10 6 G and magnetic field B = 7 x 10 10 G that 
is oriented parallel to the surface normal. The bottom panel shows the pulse fraction 
in different energy bands that is measured for IE 1207 (De Luca et al. 2004); what 
is particularly noteworthy is that the pulse fraction is larger at the spectral lines, and 
accounting for thermal effects in the model appears to be necessary to achieve these 
higher pulse fractions (see also Suleimanov et al. 2012). 

5. Discussion 

^From timing and spectral studies, CCOs appear to have relatively weak surface mag- 
netic fields (B « 10 10 — 10 11 G). The question arises as to whether CCOs are born with 
inherently weak magnetic fields (creation scenario) or they are born with much stronger 
fields but these fields were buried and are evolving (evolution scenario). The creation 
explanation is simple, but as discussed in Sec. [21 there are problems. There are also 
problems with the evolution scenario (see, e.g., Halpern & Gotthelf 2010). Nevertheless, 
evolution of magnetic fields seems natural, and there is evidence in favor of buried mag- 
netic fields in CCOs. For example, Shabaltas & Lai (2012) construct models with strong 
toroidal fields (B > 10 14 G) in the crust to explain the high pulse fraction of Kes 79 
(64 ± 2%; Halpern k Gotthelf 2010). Also Gotthelf et al. (2010) argue that a strong 
tangential field in the crust can explain the small hot spots seen on Puppis A, and this 
is confirmed qualitatively with magneto-thermal simulations by Vigano & Pons (2012). 

If CCOs have buried magnetic fields, then this sub-surface field is likely to be > 10 12 G. 
If burial is shallow, then the surface field is currently decaying. If burial is deep, then 
the surface field is growing rapidly and could lead to a rapid change in spin parameters. 
We see from Fig. [1] that CCOs reside in a relatively underpopulated region of P — P 
parameter space (see also Halpern & Gotthelf 2010; Kaspi 2010). Rapid spin evolution 
could mean that CCOs are moving out of this region quickly and joining the majority of 
the pulsar population at longer spin periods, higher P, and stronger observed magnetic 
fields. Thus magnetic field evolution could facilitate the unification of CCOs with other 
classes of neutron stars (see Kaspi 2010; Popov et al. 2010; Ho 2012). 

WCGH thanks Alexander Potekhin for discussions, appreciates the use of the computer 
facilities at KIPAC, and acknowledges support from the RAS and STFC in the UK. 
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